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ABSTRACT 

We carried out an extensive Rossi X-ray Timing Explorer campaign, in 1997 June, 
to study Circinus X-l during the active portion of its 16.55-d intensity cycle. The 
observations spanned 10 days, including 56% coverage for 7 d, and allowed us to find time 
segments which clearly demonstrate continuous evolution along the horizontal, normal, 
and flaring branches (HB/NB/FB) of a Z-source low-mass X-ray binary. These results 
confirm and extend the behavior we inferred from earlier observations. Here we study 
the continuous evolution of the Fourier power spectra and the energy spectra around the 
complete hardness-intensity track. 

A narrow quasi-periodic oscillation (QPO) peak, previously observed in the power 
spectra at 1.3-32 Hz, increases in frequency from 12 Hz to 25 Hz moving down a vertical 
extension of the horizontal branch in the hardness-intensity diagram. These horizontal 
branch QPOs (HBOs) occur near 30 Hz and fade in strength on the horizontal portion of 
the HB, while a broad peak in the power spectrum arises near 4 Hz. This peak becomes 
much more prominent along the normal branch and remains near 4 Hz (the normal 
branch QPOs, or NBOs). On the flaring branch, neither QPO is present and the power 
spectrum is dominated by very low frequency noise. We also found that each branch of 
the spectral track is associated with a specific type of evolution of the energy spectrum. 
We explored various models for the energy spectrum and parameterized the evolution 
of the spectrum in terms of a two-component model consisting of a multi-temperature 
"disk blackbody" and a higher-temperature (~2 keV) blackbody. We also show that 
an unusual line- or edge-like feature occurs at about 10 keV in energy spectra from the 
flaring branch and lower normal branch. This unusual feature is very similar to one seen 
on the FB and lower NB of the Z source GX 5—1. 



Subject headings: Stars:individual(Cir X-l) — stars:neutron — X-rays:stars 



department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 

2 Present address: Department of Physics, University of California, Santa Barbara, Santa Barbara, CA 93106; 
shirey@orion.physics.ucsb.edu 



1 



1. Introduction 

Circinus X-l exhibits a dramatic 16.55-d cycle of 
X-ray flaring which is believed to be the result of en- 
hanced mass transfer occurri ng near periastron of a 
highly eccentric bin ary orbit ( Kaluzicnski ct al. 1970] : 
Murdin et al. 1980). Although the identification of 



the orbital period seems to be secure, there is very 
little direct evidence concerning the masses of the 
two components of the binary and the other orbital 
parameters. Nonetheless, the compact component is 
thought to be a weakly magnetized neutron star on 
the basis of three type I X-ray bursts seen with EX- 
OS AT ( [Tennant et al. 1986| ). Further type I bursts 
have not been observed from Cir X-l since the EX- 
OS AT discovery, possibly because the source inten- 
sity has been higher during subsequent observations. 
No coherent pulsations, which would be expected to 
be present if the compact star is a strongly magne- 



tized neutron star, have been detected (Dower, Bradt 
k Morgan 1982| |Vaughan et al. 1994|). 



Quasi-periodic oscillations (QPOs) were seen at 
1.4 Hz, 5-20 Hz, and 100-200 Hz in EXOSAT ob- 
servations of Cir X-l in a bright state ( Tennant 1987] ; 



Pennant 1988), but other observations at lower inten- 
sity showed no such QPOs ( postcrbroek ct al. 1995 ) 



Based on these data, it was suggested that Cir X-l is 
an atoll-source low-mass X-ray binary (LMXB) that 
can uniquely reach the Eddington accretion rate and 
exhibit norm al/flaring branch QPOs ( NBOs/FBOs) 
at 5- 20 Hz flOostcrbrock ct al. 1995 ; k^an der Klii 
1994). Similar QPOs were observed in Ginga obser- 
vations of Cir X-l (|Makino 19931). 



Observations with the All-Sky Monitor (ASM) and 
Proportional Counter Array (PCA) on the Rossi X- 
ray Timing Explorer (RXTE) have shown that since 
early 1996 (from the beginning of RXTE monitor- 
ing) the baseline intensity level of Cir X-l has re- 
mained bright (-1 Crab, 2-12 keV; -1060 ^Jy at 
5.2 keV) at all phases of the 16.55-d cycle, with both 
dips and flares associated with phase zero of the cy- 
cle (Shirey e t al. 1996, 1998, hereafter Papers I & II; 
Shircy 1998). During non-flaring phases (intensity « 
1 Crab) of a cycle in 1996 March, the centroid fre- 
quency of a narrow QPO was observed to evolve be- 
tween 1.3 Hz and 12 Hz and was strongly correlated 
with the cutoff frequency of low-frequency noise and 
with the centroid frequency of a broad peak ranging 
from 20-100 Hz, at —13 times the frequency of the 
lower frequency QPO (Paper I). During portions of a 



more active cycle in 1997 February-March, a similar 
narrow QPO evolved between 6.8-32 Hz, while during 
other portions of the cycle, a broad QPO was observed 
at 4 Hz (Paper II). By correlating the timing prop- 



erties with fragmented branches in hardness-intensity 
diagrams, we identified horizontal, normal, and flar- 
ing branches. Thus, we interpreted the narrow 1.3- 
32 Hz QPOs in Cir X-l, including the 5-20 Hz QPOs 
observed with EXOSAT, as horizontal-branch oscil- 
lations (HBOs), and the broad 4 Hz QPOs as NBOs, 
thus calling into question the earlier interpretation of 
the QPOs in Cir X-l. 

In this paper we present results from a 10-day set 
of high-efficiency RXTE observations, including 56% 
coverage for 7 d starting 1.5 d before phase zero. Por- 
tions of these data show a complete Z-source track 
for Cir X-l. These results confirm our previous in- 
terpretation based on incomplete and shifted spec- 
tral tracks constructed from shorter and more widely 
spaced observations. We show how both the Fourier 
power spectra and the energy spectra evolve along the 
various spectral tracks. We explore possible models 
for the energy spectrum and parameterize the evolu- 
tion by fitting the data at various points along the 
hardness-intensity track to a model consisting of a 
multi-temperature "disk blackbody" and a higher- 
temperature isothermal blackbody. 

2. Observations 



The PCA light curves and a hardness ratio ("broad 
color") for the 1997 June observations are shown in 
Figure These data show only moderate variability 
for the first day of observations. The source entered a 
phase of significant dipping during the half day before 
phase zero (day 611. 5), based on the radio ephemeris 
( Stewart et al. 199l| ) . The hardness ratio shows that 
dramatic spectral evolution, both hardening and soft- 
ening, occurs during these dips. 

By phase zero the main dipping episode ends, and 
the transition to the flaring state begins. While the 
intensity increases by more than a factor of three in 
the lowest energy band (2-6.3 keV), the intensity be- 
tween 6.3 keV and 13 keV does not climb at all, and 
above 13 keV it decreases by a factor of about 10 
over the first 1.5 days following phase zero. This anti- 
correlation of the low and high-energy intensity dur- 
ing the transition results in a decreased hardness ratio 
after phase zero, as is observed in RXTE ASM data 
(Papers I & II; |Shirey 1998|). 
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After a relatively smooth transition toward high 
total intensity (2-18 keV) during the first day fol- 
lowing phase zero, the intensity becomes highly vari- 
able (i.e., the "active" or "flaring" state) for the re- 
maining 7 days of the observation. The intensity 
in the hard band remains much lower than before 
phase zero but shows strong variability with a peak 
value that gradually increases from about 24 to 36 



counts s PCU . Based on the results in Paper II 



and those discussed below, we can identify short-term 
variability within restricted limits as motion along 
branches in hardness-intensity diagrams and the grad- 
ual evolution of the limits as shifting of the branches. 

3. Complete Track in Hardness-intensity and 
Color-color Diagrams 

The color-color and hardness-intensity diagrams 
(CDs/HIDs) for all data in Figure @ are shown in Fig- 
ure ||. Only data from three of the five PCA detectors, 
PCUs 0, 1, and 2, were included in the diagrams, since 
PCUs 3 and 4 were not operating during portions of 
the observations. These data cover a significant por- 
tion (10 d) of an entire 16.55 d cycle. Data from 
before day 612.5 generally had soft color > 1.5 and 
broad color > 0.325, while the data after day 612.5 
generally fell below those values. 

The dips seen in Figure [l] appear as prominent but 
less dense tracks with two sharp bends in the CD (ini- 
tially toward the right of the main arc-shaped locus) 
and one sharp bend in the HID (/ < 2.3 kcounts 
s _1 PCU -1 ). We show in another paper (Shirey, 
Levine, & Bradt 1999) that tracks with these shapes 
are due to a variably absorbed bright spectral com- 
ponent plus an unobscured faint component. Brandt 
et al. (1996) used a similar model to explain the spec- 
tral changes during an intensity transition of Cir X-l 
observed with ASCA . Having identified absorption 
dip signatures, we now focus on non-dip spectral be- 
havior, which is presumably more directly related to 
the mechanisms of X-ray production. 

Most of the data in Figure || fall along a single arc- 
shaped locus in the CD and a more complicated struc- 
ture in the HID. We use spectral bands with lower en- 
ergies than those used for the 1997 February-March 
observations in Paper II; this enhances the branch 
structure in the CD. When the diagrams for the 
two observations are constructed in the same manner, 
with the same hardness ratios and only three PCUs 
rather than five, the data from these two cycles cover 



approximately the same extent in both diagrams. The 
high-efficiency coverage of the current observations 
resulted in tracks that are more complete than the 
tracks of the earlier observations. Unpublished data 
from several other cycles observed with the PCA also 
fall in similar regions of the diagrams as the data in 
Figure g. 

The detailed structure within the overall locus of 
CD/HID points is revealed in CD/HID plots of data 
divided into shorter time segments (J; 12 h). In partic- 
ular, four time segments labelled "A" , "B" , "C" , and 
"D" in Figure |l| have been selected for further study. 
The intensity from time segments A and B (days 
609.93-610.16 and 610.66-610.90 respectively) shows 
minimal variability in each energy channel, and thus 
these time segments produce small CD/HID clusters 
whose locations are indicated in Figure |2j. The source 
was more active during time segments C and D (days 
612.625-613.125 and 616.075-616.600 respectively), 
and the CD /HID tracks from these times show several 
connected branches. Enlarged views of the CD and 
HID for these time segments are shown in Figure |[ 
The full ranges of the diagrams of Figure |3 are indi- 
cated by dashed rectangular boxes in Figure 2j. Tracks 
of other time segments generally each resemble some 
portion of the entire pattern shown in Figure |, but 
often with a shifted position in the diagrams. 

The data from segment C in Figure ^ include some 
absorption dips, which result in tracks moving off the 
right side of the CD and the left side of the HID (and 
far beyond the limits of the plot in both cases). The 
timing and spectral data associated with these ab- 
sorption dips will be omitted from analysis of the HID 
branches. 

The HID patterns reveal the shape of the full 
"Z" track as previously inferred from the fragmented 
tracks in Paper II. Time segments C and D both 
show a horizontal branch, a normal branch, and a 
flaring branch which turns above the normal branch. 
In addition, segment C exhibits a long nearly-vertical 
extension on the left end of the HB, while for seg- 
ment D, there is only a small hint of an upward turn 
at the left end of the HB. The HIDs show signifi- 
cant shifts of the HB and upper NB between the time 
segments C and D, which were separated by several 
days. 

The branches in the CDs of Figure || are less well- 
separated than those in the HIDs . This was also the 
case for the diagrams in Paper II . However, the flar- 
ing branch clearly turns above the normal branch in 
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the lower left part of the current CDs, and the up- 
turned left extension of the HID horizontal branch of 
segment C is marked by an increase in the slope in 
the upper right part of the associated CD. 

The HID for segment C is similar to that derived 
from RXTE PCA observations of the Z source Cyg X- 
2 ( |3male 1998| ). The Cyg X-2 HID also shows a very 
prominent vertical extension of the HB. A similar 
upturned HB was reported in Ginga and EXOSAT 
observations of GX 5—1 ( Lewin et al. 1992] ; Kuulkers 
et il. 1994). The upturned flaring branch is similar to 
the flaring branch observed in the CD for the Z source 
G X 349+2 in recent RXTE PCA observations (^hang 
et ^,1. 199S| ). 

The HID track for time segment C was divided 
into 20 regions which were used to group data for 
further timing and spectral analysis. The 20 regions 
have been numbered as shown in Figure ^, with num- 
bers increasing from the vertical HB, through the 
NB, to the FB. We show below that region 6 does 
not adhere to the otherwise monotonic variation of 
spectral/temporal characteristics with region num- 
ber. This region may be an indication of an upward- 
shifted horizontal portion of the HB. 

Details of the temporal variability of the intensity, 
hardness ratio, and HID region numbers for time seg- 
ment C are shown in Figure |[ During this half-day 
segment, the source generally moves from lower to 
higher region number as the observation progresses. 
Thus, the time series can be divided into four sub- 
segments which predominantly correspond to each 
portion of the HID track: the vertical and horizontal 
portions of the HB, the NB, and the FB. Absorption 
dips occur in all but the flaring branch during this 
particular data set; these are easily identified by brief 
intensity dips coupled with pronounced increases in 
broad color. No region number was assigned to most 
data points from dips since the HID regions in Fig- 
ure [| were selected to avoid dips. 

The light curves of the different branches (Fig. ||) 
exhibit the following characteristics, excluding the be- 
havior associated with the dips. When the source is 
on the vertical portion of the HB, the intensity evolves 
relatively smoothly, with a slight increase in the 2- 
6.3 keV band and a decrease of almost a factor of two 
in the 13-18 keV band. On the horizontal portion 
of the HB, the source shows a substantial increase in 
soft intensity and on average shows relatively steady 
hard intensity. On the normal branch, the intensity 
is high in the soft band while decreasing and highly 



variable in the hard band. The NB/FB transition oc- 
curs at lower intensity in all bands compared to most 
of the NB. The flaring branch itself is then produced 
by high-variability "mini-flares" or bursts above the 
NB/FB apex level (region 17). 

Although the HID regions were defined such that 
obvious absorption tracks were avoided, one brief dip, 
at day 612.98, occurred from region 12 on the normal 
branch and placed a few points artificially across re- 
gions 9, 7, and 5. These points are easily identified in 
Figure ^| (bottom plot) and are thus not included in 
subsequent timing and spectral analysis. 

Likewise, the highest mini-flares on the flaring 
branch actually extend beyond region 20 and cross 
regions 8 and 9. In fact a few such points can even 
be seen above region 10 in Figure |^. The FB points 
that fell into HB regions can also be clearly identi- 
fied as points with region numbers of 8 or 9 in the 
FB portion of Figure ||. These are not included in 
subsequent timing and spectral analysis. 

4. Evolution of the Power Density Spectrum 

Fourier power density spectra (PDSs) were com- 
puted for each 16 s of time segment C. Each trans- 
form used 2 16 244-/is (2~ 12 s) time bins and covered 
the full 2-32 keV energy range. The expected Pois- 
son level, i.e., the level of white noise due to counting 
statistics, was estimated taking into account the ef- 



fects of deadtime (Morgan et al. 1997; Zhang et al 
1995; Zhang et al. 1996) and subtracted from each 



PDS; this method tends to slightly underestimate the 
actual Poisson level. For each of the 20 HID regions 
defined in Figure f|, an average PDS was calculated 
from the power spectra corresponding to points in 
that region. The average PDSs were then logarithmi- 
cally rebinned and are shown in Figure ^|. 

The general features of the power spectra are sim- 
ilar to those observed in previo us PCA observations 
(see Papers I & II). In Paper II , based on features of 
the power spectra associated with fragmented spec- 
tral tracks, we identified many of these features with 
those of a Z-source, and we discussed the properties 
of the QPOs on the horizontal and normal branches. 
In the current observations, the characteristics of the 
time variability are seen to evolve smoothly during 
a single 12-h observation (segment C). The narrow 
QPO is observed to evolve from 12 Hz in region 1 to 
25 Hz in region 7 as the HID location moves down 
the vertical extension of the HB (region 6 may be a 
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shifted version of 8). Across the horizontal portion of 
the HB (regions 8-11), the narrow QPO fades into a 
knee close to 30 Hz, while the broad QPO gradually 
rises near 4 Hz. The broad QPO is present near 4 Hz 
along the normal branch (regions 12-16). It is most 
prominently peaked in the middle of the branch and 
weak at the bottom of the branch (region 17). On the 
flaring branch (regions 18-20), no QPOs are present 
and the power spectrum shows only strong very low 
frequency noise. 

The PDS properties of the horizontal portion of 
the HB are somewhat similar to those of the upper 
normal branch of some Z sources, namely no signifi- 
cant evolution of the HBO frequency and weak NBOs. 
However, we will continue to refer to this branch as 
part of the the HB since other Z sources also show 
both vertical and horizontal portions of the HB (see 
above). 

5. Evolution of the Energy Spectrum 

The Standard2 data mode of the PCA instrument 
produces 129-channel energy spectra every 16 s. A 
parallel background file was constructed using the 
"pcabackest" program^] provided with the FTOOLS 
analysis package (version 4.0). Average pulse- height 
spectra (and background spectra) were constructed 
for each of the 20 HID regions, separately for each of 
the five PCUs. Version 2.2.1 response matrices were 
used in the analysis of these spectra. A 1% system- 
atic error estimate was added in quadrature to the es- 
timated statistical error (1 a) for each channel of the 
spectra to account for calibration uncertainties. Al- 
though the instrument response matrix is imperfectly 
known, we can safely assume that any spectral fea- 
tures that vary during the 12-hours spanned by time 
segment C are due to evolution of the source spec- 
trum. Representative spectra from the hard, bright, 
and soft extremes (regions 1,11, and 17, respectively) 
of the evolution along the HID track are shown in Fig- 
ure 0. 

The evolution of the spectrum may be studied 
by inter-comparison of ratios of pulse-height spectra 
from each region to that of a reference spectrum, from 
region 11 (Fig. ||). The spectrum is hardest in re- 
gion 1, at the top of the vertical extension of the hor- 

The background model was denned in three files provided by 
the PCA instrument team at NASA/GSFC: 

pca_bkgd_q6_e03v01 .mdl, pca_bkgd_xray_e03v02 .mdl, and 
pca_bkgd_activ_e03v03.mdl. 



izontal branch. Motion down the branch (softening, 
regions 1-7) corresponds to pivoting of the spectrum 
about ^7 keV, i.e., increasing intensity below ~7 keV 
and decreasing intensity above that energy. Motion to 
the right across the horizontal portion of the HB (re- 
gions 8-11) corresponds to continued increasing low- 
energy intensity with modest softening, but with a 
nearly constant spectrum above 12 keV. Note that in 
the hardness ratios of Figure ^ the high-energy chan- 
nel (6.3-13 keV) is dominated by photons near the 
lower bound of the interval. 

When the source moves down the normal branch 
(regions 11-17), the flux generally decreases across 
the entire 2.5-25 keV band but decreases most sig- 
nificantly at high energy from region 11 to 15 (thus 
further softening). Moving down the NB, the spec- 
trum gradually develops a dip or step above ~9 keV 
and a peak slightly above 10 keV. Motion up the FB 
(regions 18-20) is produced by increasing intensity at 
intermediate energies with a relatively constant spec- 
trum at low energy, thus hardening. The peaked fea- 
ture near 10 keV becomes more prominent moving up 
the FB and will be discussed below. 

6. Selection of Spectral Models 

Spectral forms (e.g., blackbody emission, a power 
law, etc.) for use in fitting the spectra from the 
HID regions were explored by first applying them 
to high-quality spectra from time segments A and B 
(see Fig. |l|). Variability in both of these segments 
was limited to less than 10% in all energy bands be- 
tween 2.5-18 keV. We thus constructed a single (aver- 
aged) pulse- height spectrum, known hereafter as spec- 
trum A or spectrum B, for the entire 17-19 ks of each 
segment. Errors in these spectra are dominated by 
the 1% systematics at all energies up to ~20 keV. 

The timing properties measured throughout data 
sets A and B indicate that time segment A falls on 
the vertical portion of the HB (strong narrow QPO 
at 8.4-11.5 Hz) and time segment B falls near the 
HB/NB apex (weak narrow QPO above 30 Hz and/or 
the broad 4 Hz QPO). The locations of time segments 
A and B in Figure || indicate that the branches are 
significantly shifted relative to those of segments C 
and D. 

Remillard et al. (1998) studied version 2.2.1 PCA 
response matrices using Crab nebula data and found 
that the response model was most accurate for PCUs 
0, 1, and 4 (of the five PCA detectors) and for energies 
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between 2.5 keV and 25 keV. Thus, in fitting spectra 
we only include data from PCUs 0, 1, and 4 and from 
energy channels corresponding to 2.5-25 keV. Spec- 
tra from each of these detectors are fit separately. Fit 
parameters reported are the average values for PCUs 
0, 1, and, when appropriate, 4. Errors are conserva- 
tively estimated as the entire range encompassed by 
the 90% confidence intervals from each of the detec- 
tors. PCU 4 consistently gives lower normalizations 
for fitted spectral components, so fit parameters from 
that detector are not included when computing the 



([Brandt ct al. 1996|; |Prcdchl fc Schmitt 19951) 



average normalizations and tlux values and their cr- 



rorsl Spectrum hi was not constructed lor F(JU 4 
since that detector was turned off during part of time 
segment B. 



. nterstcllar photoelectric absorption was included 



in all m odels. The absorption mo del used solar abun- 
dances ( Anders fc Ebihara 1982j ) and cross-sections 
given by Morrison & McCammon (1983). 

Several single-component models were fit to spec- 
tra A and B. Blackbody and power-law models fit 
very poorly in both cases, as did a multi-temperature 
"disk blackbody" spe ctrum (|Mitsuda et al. 198^ ; 
Makishima et al. 1986| ; model "diskbb" in XSPEC), 



with reduced x 2 (Xr) values of 22-545. A ther- 
mal bremsstrahlung model provided a better fit to 
spectrum B (x 2 = 4.0), but fit spectrum A poorly 
(Xr = 34). A relatively good fit was achieved for 
both spectra with a modified bremsstrahlung model 
(see Table [l]) which includes the effects of Compton 
scattering of bremsstrahlung photons to higher energy 
in an optically thick plasma cloud (Lamb & Sanford 
19^§ model "compLS" in XSPEC). 

A number of two-component models were also fit 
to these two spectra. A model using a disk blackbody 
and power law did not fit well (x 2 =3-5), mainly be- 
cause a single power-law slope does not adequately 
describe the spectrum at high energy. A blackbody 
with fcT> 2 keV is often included in models of the 
hard X-ray emission of LMXBs thought to contain a 
neutron star, where emission from or near the surface 
might produce high-temperature blackbody emission 
with a small effective area. Two blackbodies (kT ~ 
1.1 keV and 2.2 keV) fit moderately well (see Ta- 
ble [l]) , but required negligible interstellar absorption. 
The low absorption is inconsistent with previous mea- 
surements from ASCA and ROSAT (both sensitive 
below 2 keV where the absorption is most easily con- 
strained) which were used to estimate the interstel- 
lar column density to be N H =(1.8-2.4) x 10 22 cm -2 



Two models commonly used to fit Z-source en- 
ergy spectra are the "Western model" and the "East- 



ern model" (Hasinger et al. 1990; Asai et al. 1994). 



The Western model consists of blackbody emission 
from the hot surface of the neutron star or from 
a boundary between the accretion disk and surface, 
plus a Boltzmann-Wien component due to unsatu- 
rated Comptonization of soft photons by hot elec- 
trons (White et al. 1986; Schulz, Hasinger, & Trumper 
1989; Langmeier, Hasinger, & Trumper 1990; Bchulz 
& Wijers 1993). The Eastern model also includes 



blackbody radiation emitted from or near the surface, 
plus emission from a multi-temperature accretion disk 



ct al. 1995|) 



(Mitsuda et al. 1984; Hoshi fc Mitsuda 1991 



Hirano 



The Western model, a power law with a high- 
energy exponential cutoff plus a blackbody, fit well 
(see Table |l|) , but the best-fitting high-energy cutoff 
energy (E cut « 1.7 keV) was so low relative to the 
PC A bandpass (^2 keV) that the power law photon 
index was not well constrained. The Eastern model, 
a disk blackbody with temperatures at the inner edge 
of the disk of 1.5-1.8 keV, plus a ~2 keV blackbody, 
fit spectra A & B quite well (see Table |l| and Fig. [|) 
and gave absorption column densities roughly consis- 
tent with the ASCA and ROSAT values. Although a 
number of other two-component models also produce 
similar quality fits, the Eastern model is used below 
to provide a physically motivated parameterization of 
the spectra from the HID regions. 

The Eastern model fit to spectrum A (Fig. ^) shows 
peaked residuals at 6-7 keV, suggesting the presence 
of an emission line, probably iron Ka. Very similar 
residuals appear in most of the fits discussed above for 
both spectra A and B. Addition of a Gaussian line 
to the models does in fact improve the fits in almost 
all cases; however, the best-fitting line often has an 
extremely large Gaussian width (a > 1 keV) . The en- 
ergy resolution of the PCA is about 1 keV FWHM at 
6 keV; thus it is difficult to place reliable constraints 
on parameters such as the centroid and width of a 
narrow component. We have not included an emis- 
sion line in the fits reported in Table [l]. The pres- 
ence of an emission line near 6.4 keV is discussed in 
more detail in Shirey et al. (1999) in conjunction with 
spectra of absorption dips, which show the line more 
prominently 
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7. Fits to Spectra from 20 HID Regions 

A disk blackbody plus isothermal blackbody model 
was fit to the average spectrum for each of the 20 HID 
regions. Representative fits and residuals are shown 
in Figure The resulting fit parameters are listed in 
Tab le |2| and plotted versus HID region number in Fig- 
ure |ll| . The distance to Cir X-l has been e s timated 



to be about 6-10 kpc (Stewart et al. 1991; Goss fc 



Me jold 1977 ) , so we adopt a value of 8 kpc in con- 
verting blackbody and disk blackbody normalizations 
to radii. 

Although both spectra A and B are fit well by the 
Eastern model, the reduced chi-squared values in Ta- 
ble H indicate that none of the fits for the 20 regional 
spectra are formally acceptable. The fit results must 
therefore be regarded as an approximate description 
of the spectrum and its evolution. We emphasize that 
here, as in any case, caution is advised in drawing 
physical conclusions from the best-fit model parame- 
ters. 

The spectra along the horizontal branch (regions 
1-11) were all fit relatively well. The residuals for 
these fits (see Fig. |l^) are similar in structure to those 
for spectra A and B above. Thus they also suggest the 
presence of an emission line from iron. These spectra 
all show column densities of 1.8-2.3xl0 22 cm -2 , con- 
sistent with the ASCA and ROSAT values discussed 
above. The temperature of the ~2.0 keV blackbody 
is relatively stable on the HB (see Fig. |ll|). The tem- 
perature of the inner disk decreases from region 1 to 
region 5 (down the vertical portion of the branch) and 
then stabilizes at ~1.3 keV. The cooling is at least in 
part responsible for the pivoting of the spectrum on 
the vertical portion of the HB. 

The inner radius of the disk blackbody component, 
times a factor of order unity involving the inclination 
angle of the disk, increases from 19 to 33 km, while the 
radius of the blackbody remains between 3 and 4 km. 
These size scales are consistent with the hypothesis 
that these components arise from emission close to 
a neutron star. In this model, an increasing inner 
radius of the accretion disk is the most significant 
factor in producing the HB track; however, one should 
use caution in interpreting this as an actual physical 
radius. The inclination angle of the disk in Cir X-l is 
unknown but might be high since absorption dips are 
observed. 

From region 1 to 11, the total 2.5-25 keV flux in- 
creases monotonically, with the exception of region 6, 



from 2.89xl0~ 8 to 4.35xl0~ 8 erg cm~ 2 s _1 , corre- 
sponding to 1.2-1.8 times the Eddington luminosity 
limit for a 1.4 M Q neutron star at 8 kpc. 

Along the normal branch, the quality of the fits 
decrease from region 12 to region 17, as indicated 
by increasing Xr values (see Table ||). The absorp- 
tion column density gradually decreases by a factor 
of two, but this may be related to the decreasing fit 
quality. The inner radius and temperature of the disk 
blackbody change only slightly on the normal branch. 
In contrast, the ^2 keV blackbody begins to fade on 
the upper portion of the normal branch (regions 12- 
14), as indicated by a decreasing radius for the emis- 
sion area. The fading blackbody is illustrated in Fig- 
ure which shows the modeled incident spectra and 
both model components for several spectral fits. By 
the middle of the normal branch, the ~2 keV black- 
body has faded entirely and fits have lower % 2 values 
without it. Thus, the blackbody is omitted from the 
fits for regions 15-20. The residuals below ^6 keV 
continue to appear similar to those on the HB (see 
Fig. [n^). The peak at ~6.5 keV becomes broader and 
more complicated than on the HB, and the dip and 
peak above 8 keV become more prominent. 

On the flaring branch (regions 18-20), the fit qual- 
ity decreases further, accompanied by very low values 
for the absorption column density. A number of other 
spectral models were fit to the HID-region spectra, 
and all failed to satisfactorally fit the spectra from 
the lower portion of the HID track (region number 14 
and greater). A significant contribution to the high \ 2 
values on the flaring branch is due to the feature near 
10 keV. Addition of a Gaussian line or an absorption 
edge at 9-11 keV does improve the fits somewhat, but 
these components cannot account for all the residu- 
als near that energy. A combination of a line and 
an edge near 10 keV can adequately fit the residuals, 
but such features are difficult to justify physically at 
that energy. Even hydrogen-like iron can be ruled 
out as a possible cause due to the high energy of the 
feature. Many X-ray pulsars show cyclotron absorp- 
tion features at high energy. Inclusion of a cyclotron 
absorption component in the spectral model results 
in a fit similar in quality to that of an absorption 
edge. However, such features require magnetic fields 
of ~ 10 12 G, which would be expected to result in 
strong pulsations rather than Z or atoll behavior. 



7 



8. Discussion 



ct al. 1994), similar to our results for Cir X-l, but 



I jur spectral and timing analysis of the current ob- 



was found to be higher, 4-6 keV, in Cyg X-2 (Hasinger 



servations shows clear evidence for Z source behav- 
ior in Cir X-l. This is significant because Cir X-l 
was report ed to exhibit atoll source behavior at lower 
intensity ( |Oosterbroek et al. 19951) . Earlier RXTE 
observations, each lasting about two hours and sepa- 
rated by about two days, showed fragments of one or 
two spectral branches in hardness-intensity diagrams 
(Paper II). In the much more extensive observations 
presented in this paper, we have found longer 12-hour 
segments (time intervals C and D) which clearly ex- 
hibit all three branches of a Z source. We have demon- 
strated that these complete Z tracks shift in the HID, 
confirming the behavior we inferred from the frag- 
mented tracks of the previous observations. 

The current data also allow us to demonstrate how 
the timing properties evolve along the complete HID 
track of Cir X-l and confirm our original identifica- 
tion, in Paper H, of horizontal and normal branch 
QPOs. Fourier power spectra for different regions of 
the complete HID track show continuous evolution 
from the narrow QPO (increasing in frequency from 
12 Hz to 30 Hz in the current observations) on the 
horizontal branch, to the broad 4 Hz QPO on the nor- 
mal branch, to only very low frequency noise on the 
flaring branch. Properties of the fast timing charac- 
teristics associated with spectral branches in Cir X-l 
were discussed in Paper II . For the remainder of this 



discussion we focus on the properties of the energy 
spectrum. 

We tried fitting energy spectra of Cir X-l with 
various simple models. The spectra for time inter- 
vals A and B were well fit using the Western and 
Eastern models (see discussion below) , but no simple 
spectral form was found that fit the range of spectra 
seen during time interval C. We have not attempted 
to go beyond simple parameterized spectral models, 
e.g., by computing model spectra based on the "uni- 
fied model" of Lamb and collab orators ([iamb 1989 ; 
Psaltis, Lamb, & Miller 1995; [Psaltis ct al. 1998Q, 



which was proposed to explain the X-ray spectra and 
rapid variability of Z sources. Such sophisticated 
models may be necessary to correctly interpret the 
spectral changes in Cir X-l and other Z sources. 

The fits of spectra A and B with the Western model 
yielded a cutoff energy of ~1.7 keV for the Comp- 
tonized (Boltzmann-Wien) component. The cutoff 



ct al. 199C ). We did not use the Western model in pa- 



rameterizing evolution associated with the "Z" track 
because the cutoff energy in Cir X-l is low relative to 
the PCA bandpass, resulting in a poorly constrained 
power law index and absorption column density. 

Parameters for the Eastern model were more well- 
constrained (see Table |l|), and thus this model was 
used to parameterize the spectral variations associ- 
ated with the hardness-intensity track. In this model, 
motion along the HB is mainly associated with an 
increasing inner radius of the disk (increasing disk 
blackbody normalization) but also by a decreasing in- 
ner disk temperature. This would imply that, as the 
luminosity increases across the HB, the inner edge 
of the disk is pushed further away from the surface. 
It is not clear how this is related to the increasing 
QPO frequency, which would typically be expected 
to require a decreasing radius if the QPOs were re- 
lated to Keplerian motion at the inner edge of the 
disk, e.g., through the magnetospheric beat frequency 
model (|Alpar fc Shaham 1985|; |Lamb et al. 1985p. 



Fits of the Eastern model to energy spectra from 
the normal branch indicate that the ~2 keV black- 
body gradually fades away, leaving only the disk 
blackbody. This is similar to the result obtained when 
the Eastern model was fit to the spectrum of Cyg X- 
2, where the blackbody luminosity decreases from the 
HB to the FB (Hasinger et al. 1990). Furthermore, 
the FB of GX 5—1 is characterized by intensity dips 
which in the Eastern model can be explained by dis- 
appearance of the blackbody component, suggesting 
that accretion flow onto the neutron-star surface is 



energy in GX 5-1 was found to be 1-3 keV (Asai 



interrupted (Mitsuda et al. 1984). 

On the lower NB, a feature in the spectrum devel- 
ops above 10 keV. This feature becomes more promi- 
nent on the flaring branch. A very similar line-like 
feature at ~10 keV was reported in Ginga observa- 
tions of the Z source GX 5—1 (Asai et al. 1994). In 
GX 5—1, as in Cir X-l, the feature was present on 
the lower NB and stronger on the FB; we thus sug- 
gest that these features of the two sources may be 
of similar origin. Asai et al. showed that a peak 
near 11 keV occurs in the correlation coefficients of 
the time-series data of different energy bands versus 
the 1.7-4.0 keV band. This demonstrates a tempo- 
ral character in the narrow band at 11 keV different 
than that at adjacent energies. In turn, this gives 
assurance that the line-like feature at that energy is 
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not the result of the continuum model used to fit the 
spectrum but is intrinsic to the source. We carried 
out similar cross-correlation analysis for each of the 
20 HID regions in our study. The cross-correlation 
results, relative to the 2.5-2.9 keV band, from three 
representative regions are shown in Figure [l3|. We 
find a clear peak in the cross-correlation coefficient at 
about 11 keV in regions 19 and 20 of the FB, further 
confirming the similarity of the spectral features in 
Cir X-l and GX 5—1. Regions 16 and 17 of the lower 
NB show an abrupt drop in the cross-correlation co- 
efficient above 8-10 keV. The cross-correlation coef- 
ficients on the HB and upper NB, where the line- or 
edge-like 10 keV spectral feature is absent or weak, 
generally show much less remarkable behavior. 

Asai et al. found that that the spectral feature near 
10 keV in GX 5—1 was better fit with a Gaussian line 
rather than, e.g., an absorption edge. However, as we 
mentioned above, such a feature cannot be produced 
by even hydrogen-like iron. Asai et al. discuss several 
of the mechanisms that could possibly produce a line 
at ~10 keV. For example, emission from a heavy ele- 
ment such as Ni could produce the line, but is unlikely 
since iron should be far more abundant than the heav- 
ier elements. Alternately, a line could be blue-shifted 
due to motion in a relativistic jet or rotation in the 
accretion disk, but extreme conditions would be re- 
quired to boost the energy up to ~10 keV and one 
might expect a red-shifted line to also be observed 
in the X-ray spectrum. As mentioned above, we find 
that in Cir X-l an emission line component alone is 
insufficient to produce the observed structure of the 
feature. Both line- like and edge-like components may 
be required to explain this unusual spectral feature at 
high energy. 

Cir X-l exhibits a number of unusual features in 
its Z-source behavior in addition to the 10 keV spec- 
tral feature. Its HBOs are observed at frequencies 



as low as 1.3 Hz (Paper I), an order of magnitude 



lower than those of other Z sources. The highest fre- 
quency reached by HBOs in Cir X-l, at the HB/NB 
apex, is 30-32 Hz ( Paper II ), about a factor of two 
below the extreme HBO frequency in other Z sources. 
Power spectra of Cir X-l show a broad high-frequency 
peak, centered at 20-200 Hz, which shifts in frequency 
maintainin g an almost constant ra tio with the HBO 
frequency ( Paper I ; Tcnnant 1987 ). Cir X-l shows a 
long vertical extension of the HB. The entire HID 
track shows very large color/intensity shifts, possi- 
bly associated with the presumed eccentric 16.55-d 



binary orbit (Paper II). Atoll behavior has been re- 
ported at lower intensity ( Oostcrbroek et al. 1995| ). It 
is likely that some or many of these unusual features 
are related by some physical property of the system, 
indicating that Cir X-l may provide important con- 
straints on models of low-mass X-ray binaries. 
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at MIT and GSFC for their support. In particu- 
lar we thank E. Morgan, R. Remillard, W. Cui, and 
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work. We thank D. Psaltis for useful discussions re- 
garding Z-source spectral models. We also appreci- 
ate the detailed comments and suggestions provided 
by the referee. Support for this work was provided 
through NASA Contract NAS5-30612. 
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Fig. 1. — Light curves in three energy bands, the full 
2-18 keV band, and a hardness ratio (broad color) for 
PCA observations of Cir X-l from 1997 June 10-20, 
covering a 10-day period around phase zero (<p = 0; 
dashed line). JD 2450609.5 = 1997 June 10 h UT. 
The intensities at the beginning of these observations 
(day 610) are typical "quiescent" levels. Each point 
represents 16 s of background-subtracted data from 
PCUs 0, 1, and 2. Segments labeled A, B, C, & D 
were used for spectral studies. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v4.0. 
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Table 1 

Fit parameters for spectra A & B for four models. 
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a Errors quoted are 90% confidence limits for a single parameter (Ax 2 = 2.7). 
b Total 2.5-25 keV flux. 

c x 2 /dof, where dof = the number of spectral bins (52-54 per spectrum) minus the number of fit parameters 
(4—6). The range of values represents the fits for the several PCU detectors. 

d Blackbody radius for a distance of 8 kpc. 

°Power law with exponential cut-off above E cu t (a Boltzmann-Wien spectrum). 

f Inner radius of the accretion disk (times cos 1 / 2 0, where Q is the angle between the normal to the disk and the 
line of sight) for a distance of 8 kpc. 
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Table 2 

Fit parameters for HID regions 1-20 using a model consisting of a disk blackbody plus a 

BLACKBODY. 
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a Errors quoted arc 90% confidence limits for a single parameter (A% 2 — 2.7). 

b Inner radius of the accretion disk (times cos 1 / 2 0, where 9 is the angle between the normal to the disk and the line of sight) for a distance 
of 8 kpc. 

c Blackbody radius for a distance of 8 kpc. 
d Total 2.5-25 kcV flux. 

c x 2 /dof, where dof — the number of spectral bins (52—54 per spectrum) minus the number of fit parameters (3-5). The range of values 
represents the fits for the several PCU detectors. 
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Fig. 2. — Color-color and hardness-intensity dia- 
grams from PCA observations during 1997 June 10- 
20 (the entire period covered by Fig. |l|). Each point 
represents 16 s of background-subtracted data from 
PCUs 0, 1, and 2. Time segments A & B, defined 
in Fig. [j], produced small clusters whose locations are 
indicated in the diagrams. Time segments C & D 
produced extended tracks within the dashed boxes, 
which define the plot ranges of Fig. || 
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Fig. 3. — Color-color and hardness-intensity dia- 
grams from time segments C (left panels) and D (right 
panels) of Fig. [|. In the HID, horizontal, normal, and 
flaring branches (HB/NB/FB) have been identified. 
Each point represents 16 s of background-subtracted 
data from PCUs 0, 1, and 2. The ranges of these 
diagrams are indicated in Fig. [|by dashed boxes. 
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Fig. 4. — Hardness-intensity diagram from time seg- 
ment C. The HID track has been divided into 20 re- 
gions which were used to group data for construction 
of average power density spectra and energy spectra. 



Fig. 5. — Top: expanded light curves in 3 energy 
channels, and bottom: total 2-18 keV light curve, 
broad color, and HID regions for time segment C. The 
predominant spectral branch is identified for each por- 
tion of the data, based on HID region numbers. Ab- 
sorption dips (omitted from HID regions) are clearly 
identified by decreased intensity coupled with pro- 
nounced increases in broad-color. 
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Fig. 6. — Averaged and rebinned Fourier power den- 
sity spectra (2-32 keV) for each of the 20 regions along 
the HID track in Fig. [|. The estimated Poisson noise 
level has been subtracted from each PDS. 



Fig. 7.— Pulse-height spectra (2.5-25 keV; PCU 
only) from the hard, bright, and soft extremes, regions 
1, 11, and 17 respectively, of the HID track in Fig. || 
Labels for the region numbers are ordered vertically 
to match the relative intensities at the low and high- 
energy ends of the spectra. 
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Fig. 8. — Ratio of pulse-height spectra (PCU only) 
from selected HID regions to the pulse-height spec- 
trum of Region 11, showing the evolution of the spec- 
trum along the (a) horizontal, (b) normal, and (c) 
flaring branches. 
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Fig. 9. — Spectrum A and model (histogram) consist- 
ing of a disk blackbody and blackbody (see Table 0). 
The 1-er uncertainty, shown in the top panel as a dot- 
ted curve below the spectrum, is dominated by the 1% 
systematic error for all channels below ~20 keV. The 
residuals (x), divided by the la uncertainty, show a 
peak at 6-7 keV that may be due to an iron emission 
line. 
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Fig. 10. — Fitted pulse-height spectra and residuals 
(X) for selected HID regions (PCU only). The 1- 
a uncertainties for each spectrum, shown as dotted 
curves below the spectra, are dominated by the 1% 
systematic error for all channels below ~10 keV. Fit 
parameters are listed in Table ||. Towards higher re- 
gion number, the residuals show an unusual line- or 
edge-like feature near 10 keV. 
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Fig. 11. — Spectral fit parameters for the disk black- 
body plus blackbody model versus HID region num- 
ber. The plot is divided by vertical dotted lines 
into four parts corresponding to branches of the HID 
track. 
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Fig. 12. — Modeled incident spectra for HID regions 
1, 13, and 15, showing the two components (disk 
blackbody and blackbody) and the sum. The disk 
blackbody dominates at low energy, and the black- 
body generally dominates at high energy but has 
faded away entirely in region 15. The data shown are 
for PCU only. Fit parameters are listed in Table |^. 



Fig. 13. — Linear cross-correlation coefficients of 
time-series data (16 s bins) in each energy channel 
to that of the lowest-energy channel for HID regions 
7, 17, and 20. Towards higher region number, a step 
or peak develops in the cross-correlation coefficient 
near 10 keV and may be related to the feature seen 
at similar energy in the corresponding pulse-height 
spectra and PHA ratios. 
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